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Abstract 

Background: Previous researcli programmes tiave described muscle biocliemical traits and gene expression levels 
associated with beef tenderness. One of our results concerning the DNAJAl gene (an Hsp40) was patented. This 
study aims to confirm the relationships previously identified between two gene families (heat shock proteins and 
energy metabolism) and beef quality. 

Results: We developed an Agilent chip with specific probes for bovine muscular genes. More than 3000 genes 
involved in muscle biology or meat quality were selected from genetic, proteomic or transcriptomic studies, or 
from scientific publications. As far as possible, several probes were used for each gene (e.g. 17 probes for DNAJAl). 
RNA from Longissimus thoracis muscle samples was hybridised on the chips. Muscles samples were from four 
groups of Charolais cattle: two groups of young bulls and two groups of steers slaughtered in two different years. 
Principal component analysis, simple correlation of gene expression levels with tenderness scores, and then 
multiple regression analysis provided the means to detect the genes within two families (heat shock proteins and 
energy metabolism) which were the most associated with beef tenderness. For the 25 Charolais young bulls 
slaughtered in year 1, expression levels of DNAJAl and other genes of the HSP family were related to the initial or 
overall beef tenderness. Similarly, expression levels of genes involved in fat or energy metabolism were related with 
the initial or overall beef tenderness but in the year 1 and year 2 groups of young bulls only. Generally, the genes 
individually correlated with tenderness are not consistent across genders and years indicating the strong influence 
of rearing conditions on muscle characteristics related to beef quality. However, a group of HSP genes, which 
explained about 40% of the variability in tenderness in the group of 25 young bulls slaughtered in year 1 
(considered as the reference group), was validated in the groups of 30 Charolais young bulls slaughtered in year 2, 
and in the 21 Charolais steers slaughtered in year 1, but not in the group of 19 steers slaughtered in year 2 which 
differ from the reference group by two factors (gender and year). When the first three groups of animals were 
analysed together, this subset of genes explained a 4-fold higher proportion of the variability in tenderness than 
muscle biochemical traits. 
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(Continued from previous page) 

Conclusion: This study underlined the relevance of the GENOTEND chip to identify markers of beef quality, mainly 
by confirming previous results and by detecting other genes of the heat shock family as potential markers of beef 
quality. However, it was not always possible to extrapolate the relevance of these markers to all animal groups 
which differ by several factors (such as gender or environmental conditions of production) from the initial 
population of reference in which these markers were identified. 

Keywords: Predictors, Beef Tenderness, Array, Gene expression 



Background 

There is a need for new knowledge in order to develop 
animal farming systems that respond to the new and di- 
versified demands of producers and consumers. For pro- 
ducers, it is a question of increasing the yield of meat 
production per animal, i.e. getting heavy and lean car- 
casses at the expense of fat to ensure the highest poten- 
tial income. Consequently, genetic selection has been 
directed in favour of muscle development in order to 
produce lean carcasses in Europe. Recent genomic stud- 
ies have shown that this type of selection may change 
some muscle characteristics associated with beef quality 
[1,2]. Bovine genome sequencing has recently revolutio- 
nized the genomics research field in an unprecedented 
manner [3]. 

Although meat quality traits are moderately heritable, 
large selection programmes are very difficult to implement 
to improve beef quality due to the lack of recorded stan- 
dardized measures of beef quality [4,5]. However, consu- 
mers are looking for bovine meat of high and consistent 
quality. Like the consumers, producers are interested in 
producing high-quality meat products at the best price. 

Research done so far has highlighted the difficulties 
for control of beef sensory quality. Firstly, meat quality 
depends on complex and multifactorial biological 
mechanisms, and the muscle characteristics studied up 
to now only explain up to one fourth or one third of the 
variability in meat quality [6]. The origins of high vari- 
ability of muscle characteristics lie in the genetic poten- 
tial of the animals and to a greater extent in gene 
expression within muscles according to environmental 
factors. Thus, conditions of production (growth path, 
management, nutrition, etc) interact with the genetic po- 
tential of the animals to determine meat quality traits. 
Secondly, meat ageing plays a key role in the final ten- 
derness of beef. It is thus of paramount importance to 
increase our knowledge of all these mechanisms in order 
to better monitor meat quality components. 

Intrinsic quality attributes of beef, and especially ten- 
derness, therefore depend partly on the muscle charac- 
teristics of live animals, which in turn depend on gene 
expression. Expression levels of these genes as well as 
interaction between them can now be accessed through 
functional genomics (e.g. DNA chips and proteomic 



tools). Functional genomics is also expected to be help- 
ful for the beef industry that is looking for biological or 
molecular indicators that would identify live animals 
with desirable quality attributes, in order to direct them 
towards the most appropriate production system(s) (for 
a review, see [7-9]). 

The strategy so far has been to identify genes or pro- 
teins that are expressed differentially between animals 
displaying extreme quality attributes without any prior 
knowledge of the processes involved. Previous pro- 
grammes funded by APIS-GENE and the French Na- 
tional Agency for Research have described new genes 
associated with growth potential, beef tenderness or fla- 
vor ([1,2], for a review, see, [10]). One of these results 
concerning expression of the DNAJAl gene was patented 
[11]. The GENOTEND programme aimed to validate the 
relationship between expression of these genes and beef 
quality attributes, including tenderness and flavor. To 
achieve this objective, we designed an Agilent chip with 
specific probes for the bovine muscular genes known to 
be involved in muscle growth (including energy and pro- 
tein metabolism), carcass composition, fat metabolism 
and beef quality (including marbling). Our objective here 
was to confirm the use of this new genomic tool to pre- 
dict meat quality by looking at the previously identified 
relationships ([12]; for a review, see, [10]) between beef 
tenderness and expression of genes belonging to two 
families (heat shock proteins and energy metabolism). 

Methods 

Design of the GENOTEND chip 

More than 3000 genes involved in muscle biology or 
meat quality were selected from genetic programmes 
(for instance, [13], about 25% of the selected genes), 
proteomic studies (for instance, [14], for a review, see 
[15] ; about 14% of the selected genes), transcriptomic 
studies (for instance, [1,2,12], more than 30% of the 
selected genes), or from scientific publications (many 
genes have been previously listed by [7,9,10], about 20% 
of the selected genes). Other genes related to muscle 
biology were selected on the basis of the authors' expert 
knowledge. The list of publications from which genes 
were selected is indicated in Table 1. 
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Table 1 Main references for the selected genes in the GENOTEND chip 



Gene families Main references Proportion (%) 



Muscle hypertrophy and growth, differences between 
muscles, beef quality especially tenderness (our own 
transcriptomic studies) 


[1] Meat Sci., 82, 205-212; [2]. Meat Sci., 70, 267-277; [12], 
J. Agric. Food Chem, 55, 5229-5237; [16], EAAP 112, 371-377; 
[17]. J. Physiol. Pharmacol., 60, SuppI 2, 83-90; [18], 
J. Biochem. 133, 745-756 and unpublished results 


33 


Muscle fcetal ontogenesis, muscle hypertrophy and growth, 
beef quality (especially tendernesss) and post-mortem ageing 
(our own proteomic studies) 


[14], Proteomics 5, 490-500; [19]. Proteomics 5, 2571-2575; [20], 
Proteomics 8, 4236-4248; [21]. Animal 3, 980-1000; [22], 
Animal, 6, 867-874; [23]. 
Meat Sci., 78, 297-304 


14 


Muscle hypertrophy and growth, intramuscular fat deposition, 
beef quality especially tenderness (French, European and 
international genetic studies) 


[9]. Animal, 1, 159-173; [13]. Anim Genet, 40, 486-491; [24]. 

J. Anim. Sci., 89, 1-11; [25]. Animal, 4 303-319; [26], J. Anim. Sci., 

85, 2550-2669 and unpublished results 


25 


Muscle hypertrophy and growth, differences between muscles, 
intramuscular fat deposition, beef quality (Australian studies) 


[27], J Anim Sci 84 3239-3250; [28], BMC Dev. Biol., 7,95; [29] 
Mamm. Genome 16, 201-210; [30] Aust. J. Exp. Agr. 45, 809-820; 
[31] J. Anim. Sci. 87, 119-130. 


7 


Muscle hypertrophy and growth, differences between muscles, 
intramuscular fat deposition, beef quality and ageing 
(American and other European studies) 


[8]. Meat Sci., 74 3-15; [32]. 83, 2075-2086; [33]. 6, 2720-2731; 
[34]. 5, 1763-1769; [35]. Meat Sci., 81, 731-737; [35]. 
Meat Sci., 85, 515-518; [37]. Anim Genetics 34 438-444 


12 


Muscle biology 


Expertise of the authors 


9 



Whenever possible, several probes were used for each 
gene (e.g., 17 probes for DNAJAl). On average, four to 
five different probes were designed for each gene of 
interest. 

Animals and samples 

This study took place with two groups of 25 and 30 Cha- 
rolais young bulls and two groups of 21 and 19 Charolais 
steers (castrated at 3 months of age) slaughtered in two 
different years (years 1 and 2) from the French IVIUGENE 
programme (Genanimal-APISGENE 2005-2009). All ani- 
mals proceed from the same experimental farm 
(Bourges) and were progeny of the same seed stock sam- 
ple. The young bulls were reared in the same experimen- 
tal farm (Bourges) and slaughtered at 17 months on 
average. The steers were reared in the experimental farm 
"Le Pin au Haras" and slaughtered at 30 months on aver- 
age. All animals were slaughtered in the same conditions 
and without any electrical stimulation at the INRA ex- 
perimental slaughterhouse in compliance with current 
INRA ethical guidelines for animal welfare [22]. Muscle 
samples from the Longissimus thoracis were excised 
within 15 minutes after slaughter. Muscle samples were 
immediately frozen in liquid nitrogen and stored at -80°C 
until RNA extraction. 

In addition, steaks for each muscle were aged at 4°C 
for 14 days, and then frozen at -20°C. They were thawed 
rapidly under flowing water. On the following day they 
were grilled to a core temperature of 55°C and immedi- 
ately served to trained panellists who scored initial ten- 
derness and overall tenderness on 10-point scales: from 1 
(tough) to 10 (tender). A total of 10-12 trained panelists 
were used in each session. Panelist underwent 8-10 test 
sessions for training, and during these they evaluated 



meat from different muscles, type of animals, and cook- 
ing processes [12]. 

Hybridization studies 

After extraction, total RNA was quantified with a Nano- 
drop ND.IOOO spectrophotometer (Thermo Scientific, 
World Headquarters Location, Waltham, USA). RNA in- 
tegrity was evaluated with the 2100 bioanalyzer (Agilent 
Technologies, Massy, France) and the RNA 600 Lab 
Chip kit. The total RNA was amplified and labeled with 
Cyanine 3 using Agilent's Low RNA Input Linear Ampli- 
fication Kit, PLUS, One-Color (Agilent Technologies) 
following the detailed protocol described by Agilent. 
Briefly, 500 ng of total RNA was reverse transcribed to 
double-strand cDNA using a poly dT-T7 promoter pri- 
mer. cDNA products were then used as templates for 
in vitro transcription to generate fluorescent cRNA. La- 
beled cRNAs were finally purified using QIAGEN's 
RNeasy mini spin columns and eluted in 30 i^l of 
nuclease-free water. After amplification and labeling, 
cRNA quantity and cyanine incorporation were deter- 
mined using a Nanodrop ND.IOOO spectrophotometer 
(Thermo Scientific). For each hybridization, 600 ng of 
Cyanine 3 labeled cRNA were fragmented and hybri- 
dised at 65°C for 17 hours to an Agilent 8 x 15 K custom 
Oligo Microarray. After washing, microarrays were 
scanned using an Agilent DNA G2505B scanner. Feature 
extraction 9.1 software (Agilent Technologies) was used 
to assess fluorescent hybridization signals. 

Statistical analyses 

The total number of probes was 10257, of which 1614 
were control probes used for normalization. After data ex- 
traction with Feature Extraction (Agilent Technologies), 
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hybridization values were averaged for replicated probes. 
Generally, probes for the same gene gave similar results. 
For each array, normalization was applied with the median 
of the 1614 control probes whose average expression level 
did not significantly differ between animal groups. 

Normalization was then calculated per probe from the 
median of the probe obtained from all arrays. This 
ensured to use the same scale for all probes. Then a log2 
transformation was done. 

A Principal Component Analysis was performed with 
70 genes belonging to the heat shock protein family in 
the group of young bulls slaughtered year 1 (the refer- 
ence group) in order to identify the most correlated 
genes to initial and overall tenderness. Correlation coef- 
ficients of Pearson and of Spearman (relevant in this 
case due to the relative small number of animals per 
group) between gene expression levels and scores of ini- 
tial or overall tenderness were calculated using the Stat- 
istical software (Stat soft. World Headquarters Location, 
Tulsa, USA). Multiple regression analyses were per- 
formed using the Statistical software (Stat soft) with the 
heat shock protein genes the most correlated to tender- 
ness scores from young bulls of the reference group to 
test if these genes may predict variability of tenderness 
in the other groups of animals. Multiple regression ana- 
lyses were also performed by combining gene expression 
levels and muscle biochemical characteristics which were 
previously shown to be associated with shear force (a 
measure of beef tenderness), namely lipid content, insol- 
uble and total collagen contents and muscle fiber mean 
area [38]. The pooled relationship between muscle char- 
acteristics and gene expression with beef tenderness was 
calculated across animal groups after removing the gen- 
der and the year effects, i.e. by working on the deviation 
of each observation to the animal group muscle mean 
divided by the standard deviation for the group. 

A similar approach was conducted with 150 genes 
belonging to fat and energy metabolism. 

Results 

Muscle characteristics 

Sensory scores for the other animal groups are indicated 
Table 2. For the Charolais young bulls slaughtered year 1 
(the reference group), initial tenderness scores were on 



average 5.45 (± 1.004) while overall tenderness scores 
were on average 5.00 (± 1.115) with a coefficient of vari- 
ability of 18-20% (Table 2). Using the same experimental 
design, Guillemin et al. [38] observed that total lipid 
content, insoluble and total collagen contents as well as 
muscle fiber mean area were the only muscle biochem- 
ical traits which differ between classes of shear force (a 
measure of tenderness) for the Longissimus thoracis 
muscle, which confirmed previous results from Renand 
et al. [6]. These biochemical data were re-analysed by 
variance analysis. As known from the literature [39], we 
observed lower contents of intramuscular fat content in 
young bulls than in steers (21.2 vs. 34.3 mg/g wet tissue, 
P < 0.0001) and higher total collagen contents in young 
bulls than in steers (3.80 vs. 3.12 i^g of OH-proline / mg 
dry matter, P< 0.0001). Collagen solubility did not differ 
between genders as previously observed [39] but signifi- 
cantly differed between years (19.41 vs. 25.60 in years 1 
and 2, respectively, P < 0.04) as did mean area of muscle 
fiber (2680 vs. 3247 |im^ in years 1 and 2, respectively, 
P< 0.0001). 

Relationships between beef tenderness scores and gene 
expression levels in the reference group 

Following the Principal Component Analysis with data 
from young bulls slaughtered year 1, initial and overall 
tenderness scores were shown to be associated with the 
expression of some specific genes of the heat shock pro- 
tein (HSP) family. Thanks to correlation analyses, ex- 
pression levels of 6 genes were found to be negatively 
associated with beef tenderness (Table 3) with the high- 
est correlation coefficients (|r|>0.40). In particular, 
DNAJAl expression was negatively correlated with initial 
or overall beef tenderness (r = -0.48 to -0.59) in young 
bulls of the reference group as previously observed [12]. 
These 6 genes explain up to 49-50% of the total variabil- 
ity in tenderness in the reference group. 

With the same approach, five genes belonging to the 
energy metabolism family were also shown to be the 
most associated with beef tenderness (Table 4). We in- 
deed observed that the gene encoding the adipocyte- 
fatty acid binding protein {FABP4), expressed only in 
intramuscular adipocytes (and hence a marker of intra- 
muscular adipocyte numbers, [26]), fatty acid synthase 



Table 2 Meat sensory scores for the different groups of animals 



Initial tenderness 



Overall tenderness 



Animal group 


Number of animals 


Mean±SD 


Min-max 


Mean±SD 


Min-max 


Young bulls slaughtered year 1 


25 


5.45 ± 1 .004 


3.35-7.38 


5.00±1.115 


3.11-7.10 


Steers slaughtered year 1 


21 


5.82 ±0.91 2 


410-7.10 


5.39 ±0.949 


3.55-6.89 


Young bulls slaughtered year 2 


30 


4.64 ±0.544 


3.27-5.71 


4.54 ±0.621 


3.04-549 


Steers slaughtered year 2 


19 


454 ±0.643 


3.55-5.63 


447 ± 0.793 


3.23-6.05 


Initial tenderness and overall tenderness 


were scored within a gender x year on 10-point scales: from 1 (tough) to 10 (tender). 
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Table 3 Genes belonging to the heat shock protein family whose expression levels were the most correlated with 
initial or overall tenderness scores in the reference group (young bulls slaughtered year 1) 

Correlation with initial tenderness Correlation with overall tenderness 



Animal type Year Number of animals Person coefficient Spearman coefficient Person coefficient Spearman coefficient 



Hsp27* 


young bu Is 


year 1 


25 


-0.53 


-0.50 


-0.51 


-0.50 


DNAJA1 


young bu Is 


year 1 


25 


-0.50 


-0.59 


-048 


-0.59 


DNAJC3 


young bu Is 


year 1 


25 


-0.48 


-0.49 


-040 


-0.49 


HspHI 


young bu Is 


year 1 


25 


-0.47 


-047 


-043 


-0.47 


HspAS 


young bu Is 


year 1 


25 


-0.46 


-0.50 


-0.44 


-0.50 


HspBI 


young bu Is 


year 1 


25 


-0.46 


-0.43 


-042 


-0.42 


DNAJB5 


young bu Is 


year 1 


25 


-0.42 


-0.45 


-0.39 


-0.49 


HspA6 


young bu Is 


year 1 


25 


-0.41 


-0.40 


-042 


-0.37 


CRYAB 


young bu Is 


year 1 


25 


-040 


-0.38 


-040 


-0.41 



Among the 70 genes belonging to the heat shock protein family with probes on the GENOTEND chip, some genes were identified with the highest correlation 
coefficients with tenderness scores (they were the most associated with tenderness through the principal component analysis). Individual correlation coefficients 
(Pearson coefficient and Spearman coefficient) are indicated for each gene. Standard errors on the estimates of Pearson coefficients are between 0.18 and 0.19. 
Multiple regression analyses were performed with the combination of the first 6 genes (bold type) for young bulls slaughtered year 1. All together, these 6 genes 
explained up to 49-50% of the total variability in tenderness in the reference group. 

*Probes for Hsp27 were determined on the basis of previously published proteomic studies whereas probes for other genes were determined on the basis of 
previous transciptomic studies. 



{FASN), diacylglycerol O-acyltransferase gene {DGAT2) 
and PPAR gamma (a transcription factor regulating fatty 
acid synthesis) were positively associated with tender 
beef in young bulls slaughtered year 1. In addition, a 
gene involved in catabolism of fatty acid (hydroxyacyl- 
Coenzyme A dehydrogenize beta subunit (HADH-B, that 
catalyzes the last three steps of mitochondrial beta- 
oxidation of long chain fatty acids) was associated with 
the most tender beef in young bulls slaughtered year 1. 
All together, these 5 genes explain up to 66-71% of the 
total variability in tenderness in the reference group. 

Cross-validation of the identified tenderness markers 

Five groups of probes (corresponding to four genes) belong- 
ing to the HSP family (Hsp27 encoded by the HspBl gene, 
identified from proteomic studies with the accession num- 
ber Gene Bank ID: GI:71037405, HspBl Gene Bank ID: 
NM_013560.2, DNAJAl Gene Bank ID: NM_001015637.1, 



HspHl Gene Bank ID: NM_001075302.1, and HspAS Gene 
Bank ID: NM_174345.3) were found to be significantly 
expressed in all animals across genders and years of slaugh- 
ter (which means expression levels differ from zero). In 
contrast, only two genes related to fat metabolism (FASN 
and HADHB) were found to be significantiy expressed in all 
animals and were not validated for all animal groups. 
Therefore, we attempted to validate only the five groups of 
probes belonging to the HSP family in all animal groups. 

Multiple regression analysis showed that the five 
selected groups of probes belonging to the HSP family 
explain 39 and 44% of initial and overall tenderness re- 
spectively in the reference group (young bulls slaugh- 
tered year 1) in which they were identified (Table 5). 
Furthermore, they also explained from 37 to 49% of the 
variability in initial and overall tenderness in young bulls 
slaughtered year 2 and in steers slaughtered year 1 
(P < 0.08), but not in steers slaughtered year 2 for which 



Table 4 Genes belonging to the energy and fat metabolism family whose expression levels were the most correlated 
with initial or overall tenderness scores in the reference group (young bulls slaughtered year 1) 





Animal type 


Year 


Number of animals 


Correlation with initial tenderness 


Correlation with overall tenderness 


Person coefficient 


Spearman coefficient 


Person coefficient 


Spearman coefficient 


FABP4 


young bu Is 


year 1 


25 


0.42 


0.31 


041 


0.39 


PPARG 


young bu Is 


year 1 


25 


0.59 


048 


0.52 


0.46 


DGAT2 


young bu Is 


year 1 


25 


0.61 


0.51 


0.57 


0.56 


FASN 


young bu Is 


year 1 


25 


0.61 


0.56 


0.59 


0.52 


HADHB 


young bu Is 


year 1 


25 


0.44 


049 


047 


0.46 



Among the 150 genes belonging to the energy and fat metabolism family with probes on the GENOTEND chip, five genes were identified with the highest 
correlation coefficients with tenderness scores (they were the most associated with tenderness through the principal component analysis). Individual correlation 
coefficients (Pearson coefficient and Spearman coefficient} are indicated for each gene. Standard errors on the estimates of Pearson coefficients are between 0.17 
and 0.1 9. Multiple regression analyses were performed with the combination of these 5 genes for young bulls slaughtered year 1 . All together, these 5 genes 
explained up to 66-71% of the total variability in tenderness in the reference group. 
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Table 5 Validation in the different animal groups of the five sets of probes belonging to the heat shock protein family 
to predict variability in tenderness 









Initial tenderness 


Overall tenderness 


Animal group 


Number of animals 






P value 




P value 


Young bulls slaughtered year 1 


25 


HSP gene expression levels 


0.39 


0.07 


0.44 


0.04 






Biochemical traits 


0.08 


0.75 


0.08 


0.78 


Steers slaughtered year 1 


21 


HSP gene expression levels 


0.45 


0.08 


049 


0.05 






Biochemical traits 


0.42 


0.06 


041 


0.06 


Young bulls slaughtered year 2 


30 


HSP gene expression levels 


0.37 


0.04 


040 


0.03 






Biochemical traits 


0.04 


0.88 


0.08 


0.69 


Steers slaughtered year 2 


19 


HSP gene expression levels 


0.20 


0.68 


0.23 


0.58 






Biochemical traits 


0.19 


0.52 


0.21 


0.49 



Multiple regression analyses were performed with the combination of the 5 groups of probes [DNAJAt, HspHl, HspA8, plus two groups of probes for Hsp67-Hsp27) 
to predict beef initial or overall tenderness for young bulls or steers slaughtered year 1 or year 2. Multiple regression analyses were also performed with the 
combination of four muscle biochemical traits {Intramuscular fat content, mean fiber area, total collagen content and collagen solubility) previously shown [38] to 
be associated with beef quality. The proportion of the variability in tenderness which is explained by the combination of the 5 groups of probes is indicated (R^) 
as well as the associated P value. 



only 20-23% of the variability in tenderness can be 
explained (not significant, Table 5). Regression analysis 
was found to be very similar in three animal groups 
(young bulls slaughtered years 1 and 2, steers slaugh- 
tered year 1). This result validates the relevance of this 
set of genes identified in one group of Charolais young 
bulls to predict beef tenderness at least in two other 
groups which differ from the reference population by 
only one factor (animal gender or environmental condi- 
tions reflected by different years of slaughter). However, 
when animals differed from the reference population 
both by gender (castration vs. intact males) and environ- 
mental conditions (slaughter year 2 instead of year 1) 
such as in the case of steers slaughtered year 2, the 
selected genes cannot be used any more to predict beef 
tenderness at a significant level. 

With the same approach, we observed that the four 
muscle biochemical traits (intramuscular fat content, 
total and insoluble collagen contents and muscle fiber 
area), which were previously identified to differ between 



shear force classes in the Longissimus thoracis of the 
same animals [38], explained only 8% of the variability in 
tenderness in young bulls slaughtered year 1 (not signifi- 
cant). As well, they explained only 4-8% and 19-21% in 
young bulls and steers respectively slaughter year 2 (not 
significant). However, they could explain 42% of the vari- 
ability in tenderness in steers slaughtered year 1 
(P = 0.06, Table 5). 

When the three animal groups in which the HSP genes 
were validated were analysed together after removing the 
gender and the year effects (namely a total of 76 animals: 
young bulls of years 1 and 2, steers of year 1, Table 6), 
this set of genes explained 20% of the variability in the 
initial and overall tenderness (P = 0.008), whereas the 
biochemical muscle traits known to be associated with 
tenderness (namely intramuscular fat content, mean fiber 
area, total collagen content and collagen solubility 
according to [38]) explained only 5-6% of the variability 
in the initial and overall tenderness (not significant. 
Table 6). 



Table 6 Validation across genders and environmental conditions of the set of five groups of probes belonging to the 
heat shock protein family to predict variability in tenderness 







Initial tenderness 


Overall tenderness 


Predictors of tenderness 


Number of animals 


R^ 


P value 


R^ 


P value 


Expression levels of H5pfi/fHsp27), DNAJAl, HspHl, HspA8 


76 


0.20 


0.008 


0.20 


0.008 


ntramuscular fat content, mean fiber area, 
total collagen content and collagen solubility 


76 


0.06 


0.39 


0.05 


0.49 


ntramuscular fat content, mean fiber area, total collagen 
content and collagen solubility plus expression levels 
of HspBl (Hsp27), DNAJAl HspHl, HspA8 


76 


0.27 


0.07 


0.25 


0.02 



Multiple regression analyses were performed with the combination of the 5 groups of probes {Hsp27, DNAJAl, HspHl, HspA8, plus two groups of probes for HspBl- 
Hsp27) to predict beef initial or overall tenderness for young bulls slaughtered year 1 (the reference group) or year 2 or steers slaughtered year 2. The four muscle 
biochemical traits (Intramuscular fat content, mean fiber area, total collagen content and collagen solubility) previously shown to be associated with beef quality 
[38] were analysed with a similar approach alone or in association with the 5 groups of probes {DNAJAl, HspHl, HspA8, plus two groups of probes for HspBl- 
Hsp27). The pooled relationship between muscle characteristics and gene expression with beef tenderness was calculated across animal groups after removing 
the gender and the year effects. I.e. by working on the deviation of each observation to the animal group muscle mean divided by the standard deviation for the 
group. The proportion of the variability in tenderness which Is explained in the different situations is ind icated (R^). 
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Discussion 

Of all the beef quality traits, tenderness is considered to 
be the most important. It appears that among the main 
determinants of tenderness are collagen characteristics 
and the muscle fiber type [6] and, in association with 
this, the extent of proteolysis during post-mortem ageing 
and also the amount of intramuscular fat which facili- 
tates mastication of beef Indeed, with the same animal 
groups, intramuscular fat content and collagen solubility 
were shown to be positively associated with tenderness, 
whereas muscle fiber cross-sectional area were shown to 
be negatively associated with tenderness [38]. However, 
the contribution of collagen characteristics to tenderness 
may depend on the cooking temperature (55°C in our 
work, which is lower than in other studies) since heat- 
induced changes of intramuscular connective tissue were 
observed [40]. Among the proteolytic systems which 
control post-mortem proteolysis [41], it has been re- 
cently proposed that the apoptotic pathway may be of 
great importance, since it is the first event occurring 
after the animal's death. Our results confirm these 
observations using a new DNA chip. The aim of this 
manuscript was hence to validate this tool by confirming 
two well-accepted or recently discovered biological 
determinants of beef tenderness: the relationships of 
beef tenderness scores (i) with expression levels of genes 
corresponding to heat shock proteins on one hand and 
(ii) with expression levels of genes related to fat and en- 
ergy metabolism in muscle fibers on the other hand. 

Relationships between beef tenderness scores and 
expression levels of genes encoding heat shock proteins 

After slaughter of the animals, there is generalized cell 
death in all organs and tissues due to the lack of oxygen 
and hence anoxia of all cells. This stage precedes rigor- 
mortis and subsequently meat ageing during which pro- 
tein degradation occurs thanks to the activity of different 
proteolytic systems (for a review, see [42]). More and 
more evidence is now described in the literature indicat- 
ing that the muscle death process at the beginning of 
conversion of muscle into meat may play an important 
role in regulating meat quality [41]. This cell death 
process, often referred to as apoptosis, is regulated by 
many proteolytic enzymes including caspases, the activ- 
ity of these enzymes being itself regulated by many 
factors. 

Cell death induces considerable cell stress and, in re- 
sponse to stress, cells rapidly produce heat shock pro- 
teins (HSP) that play a universal role in maintaining 
cellular homeostasis. HSP help in maintaining the integ- 
rity of the cell and have an anti-apoptotic activity [43]. 
For instance, proteomic studies showed that Hsp27 in 
fresh Blonde dAquitaine muscle and levels of Hsp27 
fragments in aged meat explained up to 91% of variation 



in sensory scores [23]. It was hypothesized that higher 
levels of Hsp27 associated with limited aggregation of 
muscle proteins could facilitate the action of proteolytic 
enzymes during meat ageing. Therefore, the positive cor- 
relation between tenderness scores and levels of specific 
Hsp27 protein isoforms or fragments in Blonde dAqui- 
taine [23] contrasts with the negative correlation we 
observed in this work between Hsp27 mRNA levels and 
tenderness scores. This may be explained by the differ- 
ent sources of variability between the two studies: ageing 
time [23] or animal variability in this study. However, 
our results fit well with the positive correlation of Hsp27 
protein level and shear force value in Korean cattle [44]. 
Indeed, a recent study with French breeds confirms that 
correlation of Hsp27 level may be positive or negative 
depending on the breed [45]. Other authors have 
described the complex effects of meat pH upon the sub 
cellular distribution of muscle HSP during ageing of beef 
[46] . In other experiments, HspBl (encoding Hsp27) and 
its regulator genes were shown to be negatively corre- 
lated with intramuscular fat content while being asso- 
ciated with a high shear force [47] and hence low 
tenderness as in our work. 

Our results with the GENOTEND chip also confirm 
the negative correlation between DNAJAl expression 
and the tenderness score [12], but for young bulls and 
steers of year 1 only, which confirms our hypothesis that 
some markers of beef quality are highly dependent on 
rearing practices and environmental conditions. DNAJAl 
encodes a member of the large 40 kDa heat shock pro- 
tein family (Hsp40). This protein is a co-chaperone of 
the 70 kDa heat shock protein (Hsp70) and is believed 
to play a role in protein folding and mitochondrial pro- 
tein import. The DNAJAl IHs^lO complex also directly 
inhibits programmed cell death (or apoptosis), which 
supports the hypothesis of Ouali et al. [41] that apop- 
tosis is important for the tenderization of beef during 
ageing. Depending on the animal group (young bulls or 
steers slaughtered year 1), DNAJB9, DNAJC3 or 
DNJACIO, which are also members of the Hsp40 family, 
were negatively correlated to tenderness scores (data not 
shown) just like DNAJAl, which supports the role of 
members of the Hsp40 family in tenderness. In humans, 
the DNAJ family has over 40 DNAJ members [48] and 
can be subdivided into three subfamilies: DNAJA pro- 
teins, DNAJB proteins and DNAJC proteins with differ- 
ent structural features. Our results also support the role 
of some HspA (Hsp70), and some HspB (small HSPs) in- 
cluding Hsp27, in addition to that of DNAJAl (Hsp40). 
For instance, CRYAB, which encodes aB-crystallin and 
shares homology with Hsp27, was also negatively corre- 
lated with tenderness in steers. However, the precise 
genes of this family which are individually negatively 
correlated to tenderness scores depends on the animal 
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group and hence on its characteristics (gender, breeding 
management, etc). 

Nevertheless, a major finding of our study is that a 
subset of genes belonging to the heat shock protein fam- 
ily may be suitable to better explain than muscle bio- 
chemical characteristics the variability of tenderness 
scores for both two groups of Charolais animals of the 
same gender or reared in similar environmental condi- 
tions. However, when the animals differ from the refer- 
ence population by more than one factor (gender, or 
environmental conditions, etc), the expression levels of 
these genes cannot explain anymore the variability of 
tenderness scores. Further work is thus needed to extent 
the list of gene markers to have a more generic predic- 
tion of tenderness in other conditions of rearing and/or 
slaughtering and, if possible, in other breeds. 

Relationships between beef tenderness scores and 
expression levels of genes related to muscle metabolism 

As previously demonstrated recently by Zhao et al. [49], 
we observed that some genes known to be expressed 
only in adipocytes {FABP4), or involved in lipid depos- 
ition {PPARG, DGAT2, FASN) were characterized by ex- 
pression levels positively associated with tenderness in 
young bulls slaughtered year 1. This strongly suggests 
that intramuscular fat level contributes to tenderness in 
agreement with previous research. Indeed, whereas it is 
well-known that intramuscular fat level affects mainly 
juiciness and flavor of beef, it also affects indirectly ten- 
derness [50]. In addition, with the same samples of this 
study, Guillemin et al [38] observed that intramuscular 
fat content was 45% higher in the most tender meat 
samples (determined from the shearforce measurement) 
compared to the least tender meat samples. Intramuscu- 
lar fat is mainly stored within intramuscular adipocytes 
embedded in a connective tissue matrix in close proxim- 
ity to a blood capillary. So, in fat-rich beef, intramuscular 
fat deposited between fibers fascicules may disrupt the 
structure of muscle connective tissue, thus contributing 
to increase meat tenderness [51]. This is important in 
beef meat, in which the abundance and solubility of col- 
lagen (the main component of connective tissue) are a 
source of variation in tenderness especially in young 
bulls. Generally, intramuscular fat has a diluting effect 
on connective tissue (known to increase toughness) and 
consequently a tenderizing effect on meat. Meat from 
animals with a high intramuscular fat level often has a 
low shear force [47], and hence high tenderness. There- 
fore, it is not surprising to observe that FABP4 expres- 
sion level, known to be a marker of intramuscular 
adipocyte number [26], was positively correlated with 
tenderness. Similarly, the more fatty acids are stored in 
these adipocytes, the bigger the adipocytes are, which 
favors tenderness by diluting connective tissue. Storage 



of fatty acids is controlled by two mechanisms. The first 
is due to the action of LPL which hydrolyses circulating 
triglycerides and hence favors uptake of resulting fatty 
acids. The second is associated with i) the activity of 
FASN which is a key enzyme involved in the de novo 
fatty acid synthesis and more precisely the elongation of 
fatty acids and (ii) the activity of DGAT which catalyses 
the formation of triglycerides from diacylglycerol and 
Acyl-CoA. 

It is also well-known that muscles composed of fast 
fibers which low intramuscular fat content are more sus- 
ceptible to early post-mortem proteolytic degradation 
than muscles mainly composed of slow fibers. Fast fibers 
tend to have higher levels of stored glycogen (which 
favors pH decline during ageing) and higher amounts of 
Ca2 + -activated myosin ATPase (which favors the speed 
of ageing). Therefore, after a short period of ageing, we 
expect fast muscles to generate tender beef compared to 
slow muscles due to a faster rate of ageing (for a review, 
see [52]). However, this may be not true after 14 days of 
ageing when a significant part of differences in ageing 
speed have disappeared. This may explain why the pro- 
portion of slow fibers has been reported to improve ten- 
derness in cattle, especially for Longissimus muscle 
although this relationship may be not true for other 
muscle types or not always valid with data across muscle 
types (for a review, see [53]). However, the observation 
that oxidative fibers are associated with tender beef has 
been confirmed in our study by a higher expression of 
the HADHB known to be mainly expressed in oxidative 
fibers and also in tender beef in our case. 

In contrast to the situation with HSP genes, we found 
no subset of genes related to fat metabolism or muscle 
fiber metabolism which were significant across environ- 
mental conditions and genders to predict beef tender- 
ness. In some ways, this observation confirms results 
recently obtained by de Jager et al. [54] who observed 
differences in expression of the lipid storage genes be- 
tween tough and tender genotypes, but also observed 
large differences in the expression of genes involved in 
fat and energy metabolism across two different experi- 
mental sites in Australia independently of genotype. A 
likely explanation is that these gene expression changes 
are reflecting differences in metabolic status of the ani- 
mals such as turnover rates of nutrients (fatty-acid pre- 
cursors such as glucose, triglycerides and fatty acids) or 
mitochondrial activity, especially considering nutritional 
regimens, feed density, ambient temperature, slaughter- 
ing conditions, etc between the different experimental 
sites or the different years of slaughtering. 

Conclusions 

The GENOTEND chip was designed and validated by 
confirming known or recently discovered relationships 
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between tenderness scores and expression levels of genes 
encoding heat shock proteins or related to muscle fiber 
metabolism. Furthermore, this transcriptomic analysis, 
focused on some muscle/meat quality markers, provided 
the means for analysis of these genes in different animal 
types (young bulls or steers) reared in different environ- 
mental conditions (corresponding to two years of 
slaughtering). We observed that numerous markers of 
beef tenderness can be identified but they are often spe- 
cific to an animal type (steer or young bull) or to envir- 
onmental conditions. However, these results confirm the 
idea that heat shock proteins (especially of the Hsp40 
family) on the one hand [12,41] or metabolic enzymes 
on the other hand ([26]) may be potential good markers 
of beef tenderness. More specifically, a subset of genes 
belonging to the HSP family is able to explain a signifi- 
cant proportion of the variability of tenderness across 
genders and environmental conditions in three animal 
groups over four unlike some muscle biochemical char- 
acteristics yet known to be negatively or positively asso- 
ciated to shear force of beef 
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